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Analysis of a chemotaxis system modeling ant 

foraging 

Ricardo Alonso^, Paulo Amorim^ and Thierry Goudon^ 


Abstract 

In this paper we analyze a system of PDEs recently introduced in [P. 
Amorim, Modeling ant foraging: a chemotaxis approach with pheromones 
and trail formation], in order to describe the dynamics of ant foraging. 

The system is made of convection-diffusion-reaction equations, and the 
coupling is driven by chemotaxis mechanisms. We establish the well- 
posedness for the model, and investigate the regularity issue for a large 
class of integrable data. Our main focus is on the (physically relevant) 
two-dimensional case with boundary conditions, where we prove that the 
solutions remain bounded for all times. The proof involves a series of fine 
a priori estimates in Lebesgue spaces. 

Keywords. Ant foraging, Chemotaxis, Animal movement. Reaction-diffusion equa¬ 
tions. 
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1 Introduction 

One of the most interesting phenomena arising in the collective behavior of ants 
is the formation of trails. Indeed, while each individual ant has a very limited 
cognitive ability, the population as a whole is capable of complex, organized col¬ 
lective behavior, such as brood rearing, waste management or fungus gardening. 
Even more striking is the fact that many of these activities, especially trail for¬ 
mation (which occurs during foraging, migration, or aggression), are essentially 
leaderless and yet highly organized. 

Many tools have evolved in ant societies to allow for this sort of complex, so- 
called emergent behavior. One of the most important is the use of pheromones 
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as a means of communication between individuals. Pheromones are chemical 
compounds secreted by ants (and many other species of animals) which are used 
to convey information between individuals and to signal different states, such as 
presence of food or alarm. Each pheromone triggers a corresponding behavior 
in individuals: when an ant detects alarm pheromone, it becomes itself alarmed 
and secretes more pheromone, leading to a chain reaction among the population, 
whose effect is to elicit an apparently organized defense response. 

We are interested here in the trail-following behavior of ants, which is trig¬ 
gered in part by trails of pheromones. When an individual ant, foraging at 
random, encounters a food source, it typically travels back to the nest leaving 
on the substrate a trail of pheromone. When other ants encounter the trail, they 
follow it in the direction of the food, and upon finding the food, they head back 
to the nest and deposit more pheromone. Thus, as long as food is available, 
the trail will be reinforced and the food will be removed efficiently. Conversely, 
when the food is exhausted, the evaporation and diffusion of the pheromone 
quickly erase the trail, when it is no longer being reinforced. 

The dynamics of ant foraging behavior has recently come under increased 
interest from mathematicians trying to find suitable frameworks in which to 
analyze this behavior. We refer the reader to the recent works dull [50], 
among others. In d (see also the independent work d), a system of PDEs, see 
(ISPDl) below, is introduced in order to describe the dynamics of ant foraging. 
Roughly speaking, the population splits into two parts: the ants searching for 
food, and the ants going back to the nest, and the pheromone production is 
interpreted as a chemotactic mechanism that drives the population to privileged 
directions. The discussion on the modeling issues in d is complemented by a 
set of convincing numerical simulations that illustrate the ability of the model 
to reproduce relevant behaviors of ant populations. Here, we aim at analysing 
the mathematical properties of these equations. 

More generally, finding relevant models able to reproduce the formation of 
the space-time heterogeneous patterns observed in life sciences is becoming a 
very active field, particularly motivated by the landmark contributions of T. Vic- 
sek et al [32, and F. Cucker and S. Smale duni about the formation of flocks 
in large populations of birds or fish. The key feature relies on the transmission 
by the individuals of the information contained in their close environment, so 
that the whole population organizes according to remarkable patterns. In this 
vein, various models have been introduced, which have led to original problems 
for mathematical analysis and fascinating numerical simulations that reproduce 
certain features of natural phenomena; we refer the reader to the surveys [241I33| 
for an overview on the subject. Here, following di are interested in con¬ 
tinuum models, where populations are described by their local concentrations. 
The interaction between the individuals can be though of through a certain po¬ 
tential, which is defined self-consisteiitly, thus depending on the variations of 
the concentrations. This is reminiscent to the theory of chemotaxis, which dates 
back to C. Patlak |3S] and E. Keller-L. Segel [331 [35] to model the behavior of 
certain bacteria and slime molds, which are attracted to chemical substances 
that they themselves secrete. In particular, the possible formation in finite time 
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of singularities in the solutions of the Keller-Segel equations, the concentration 
into Dirac masses corresponding to the aggregation of the population into a 
single location, has motivated a huge amount of mathematical works, see for in¬ 
stance nsumiini- By now, the mathematical theory of the Keller-Segel system 
is well established, and we refer the reader to the surveys HZIIIHIIIII for further 
information and references. To describe interaction mechanisms between liv¬ 
ing organisms by chemotactic principles has been successfully adapted to many 
different situations, see e. g. Em [II muni 123 [SB to mention a few. 

To be more specific, in this paper, we study the basic mathematical theory 
for noiinegative solutions (t, x) i—>■ {u, w,p, c)(t, x) of the following model for ant 
foraging 

dtu — Au + V- (u yVp) = —uc + wN, 

dtw — DyjAw + V- {w Vu) = uc — wN, 

< (SPD) 

dtp — DpAp = Pw — 5p, 

dtc = —u c. 

In this system of PDEs, the unknowns are 

• the density of foraging ants u, 

• the density of returning ants w: it describes the ants which have found 
food and are returning to the nest, 

• the concentration of the pheromone p, 

• the distribution of the food c. 

These nonnegative quantities depend on time {t > 0) and space (a; G D C M") 
variables. The data of the problem are 

• the site of the nest, embodied into the function x i—>■ N{x), 

• a function x i—> Pix) that describes the pheromone deposition as returning 
ants approach the nest; typically this function decreases as the distance 
to the nest decreases, 

• a nest-bound vector field x i—>■ Vv{x) representing the speed of the ants 
when returning to the nest (it might contain information on the topogra¬ 
phy, obstacles...), 

• Diffusion, sensitivity and evaporation coefficients Dp, x, and <5, which 
are all positive constants. 

The system (ISPDI) will be addressed in the sequel as the Slow Pheromone 
Diffusion model as the pheromone diffusion time scale is comparable to that 
of the dynamics of the ant foraging. We refer the reader to [2] for details on the 
biological motivation for system (ISPDl) . In what follows, we will take P = 1 and 
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Dp = Du, = X = 1 for the sake of simplicity. We will study the system (ISPDl) 
on the whole Euclidean spatial domain with a few comments applying to 
dimension n = 3. 

We will also work with the simplified situation where the pheromone diffusion 
time scale is small compared to the dispersal of the ants. In addition, we suppose 
that there exists a very abundant, or renewable, food source 0 < c := c{x), so 
that we can assume it is a given function of space. These simplifying assumptions 
lead us to the following reduced system, with unknowns (t,x) i—>■ {u,w,p){t,x) 

' dtu — Alt + V- (u Vp) = —uc + wN, 

< dtw — Aw + V- {w Vv) = uc — wN, (FPD) 


— Ap = w — Sp. 


It will be referred to as the Fast Pheromone Diffusion model. This system 
will be considered in a domain fl C having a smooth boundary d^l. In order 
to conserve the total mass of ants, we impose the following zero-flux boundary 
conditions 


Vm • n 


an 


= (wVa — Vw) • n 


an 


= Vp • n 


an 


= 0 , 


( 1 . 1 ) 


where n stands for the outward unit normal vector to dH,. While this is not 
crucial for most of the analysis, we can bear in mind the fact that physically 
relevant velocity fields satisfy 


Va-n|g^<0 (1.2) 

since Vu is pointing towards the nest. Assumption m will play a role when 
proving uniform propagation of the L°“-norm for solutions, namely, estimate 
(I3.23P and later, for the local existence of classical solutions. Therefore, ants do 
not escape the domain Q. The initial data 


u 


= Uo, 

t=0 


w 


= Wo 

t=0 


(1.3) 


will be assumed, naturally, as nonnegative integrable functions. As a conse¬ 
quence, we have 


/ {u + w){t,x) dx = / {uo + Wo){x) dx. 

Jn Jn 


(1.4) 


We shall prove the global well-posedness in dimension n = 2 for both models 
(ISPDl) and (IFPDI) . In fact, weak solutions of such systems are bounded in 
(0,T] X fl for any T > 0. This is in contrast with the situation known for the 
usual Keller-Segel system 


dtp -I- V-(pV$) - Ap = 0, 
- A$ = p. 


(1.5) 
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It is indeed well-known that, if the initial mass f p(0,x)dx exceeds a certain 
threshold (which depends on the domain and the space dimension n > 2), then 
the solution of (II3D blows up: ||/9(f,a:)||oo —>■ oo as t T* < oo, typically 
exhibiting a concentration to a certain location. Roughly speaking, this effect 
can be seen as a result of the competing effects between diffusion, and the 
explosive behavior of the ODE y = (which is obtained neglecting diffusion 
and looking at the solution p evaluated along the characteristics associated to 
the field Vd)). At first sight, the structure of the systems (ISPDI) and (IFPDI) is 
quite close to that of the Keller-Segel system and one may wonder whether 
of not such threshold phenomena occur to produce the blow up of the solutions. 
However, we shall see that the introduction of an additional population (the 
returning ant population w), which is itself subjected to a regularizing parabolic 
equation, prevents the blow up formation. 

The paper is organized as follows. In Section [5] we set up the assumptions 
necessary to the analysis and we state the main results. Section [3] is devoted to 
the model (iFPDll and it is divided in four subsections for clarity. The first and 
second subsections present the core of the analysis where increasing integrability 
of the solutions is shown by means of a priori estimates. The initial step is to 
prove instantaneous generation of L'>'-integrability. To be more specific, we 
prove that initial data will lead to solutions lying in L'>'(D)) for 

any 7 S (l,c») and any positive times 0 < < T <00. The bounds given 

for the L'''-norm of the solution will depend on the structure of the system, that 
is, on the data N, c, v and more importantly, on the conserved quantity, that 
is the initial mass of the total ant population mo = J Uodx + f Wo da;. Such 
bounds are independent of the time existence interval, a fact that can be used 
to justify that the solutions are globally defined. Furthermore, it is possible 
to use such a result to prove that solutions are actually uniformly bounded for 
any positive time. The proof of this fact relies on the De Giorgi energy level 
set method, m- The third subsection is concerned with the global existence 
of classical solutions for the problem. Given smooth initial data, local in time 
existence of solutions can be justified by using the classical Banach fixed-point 
theorem on suitable metric spaces. Next, the global a priori estimates which are 
valid for classical solutions allow us to justify that the lifespan of these solutions 
is actually infinite. We finish Section |3] by presenting a global well-posedness 
theorem for weak solutions of (TfpdI) with initial data in fl L'*' with 7 > 2. 
The proof uses approximation by the classical solutions just found. Finally, in 
Section |4] we discuss the well-posedness of the system (ISPDil . Our approach 
uses elementary properties of the heat kernel which is, somehow, a different 
approach than the one used for the analysis of (iFPDil . 

2 Notations, hypotheses and main results 

The main results of this paper are concerned with the well-posedness of weak 
solutions of the system (iFPDil . 
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Definition 2.1. We say that the triple {u, w,p) is a weak solution of the system 
(iFPDl) if it satisfies: 

(i) {u,w) G L2(0,T;iJi(O)) and {dtu,dtw) G {Q,T-. 

(ii) Equations (iFPDl) are solved in the sense that, for any test function ^ G 
C°°([0, oo) X n), compactly supported in [0,T) x SI, we have 


/7 

/o Jq 


/7 


[ —udt( + {Vu — uWp) ■'S/C){t,x) dx dt — / Uo{a::)C(0, x) dx 

Jn 

= / {—uc + wN)({t,x)dxdt, 

Jo Jn 

{ —wdtC + {'S/w — wWv) ■'S/(){t,x) dx dt — / ■u;o(a;)^( 0 , a:) dx 

Jn 


[uc — wN)((t, x) dx dt, 


>0 Jn 


/ Vp-VCdx= / {w — Sp)Cdx. 
Jn Jn 


Main hypotheses. The following assumptions will be used throughout this 
Section: 

(i) The initial data {uo,Wo) is nonnegative with finite mass: 

/ Uo(a:) + / Wo{x) dx = nio < oo . ( 2 - 1 ) 

Jn Jn 

(ii) The parameters of the system (IFPDl) are such that 

{N, c, Vv, Av) G (H) 

The results can be extended under weaker assumptions on v at the price 
of more involved technicalities in the estimates. As a convention, let us 
agree here that when a constant is referred to depend on (El it means 
that this constant depends on the L°°-norms of N, c, Vv and Av. 

(Hi) The domain il is of class C^. 

In what follows we will use the shorthand notation 7 ^ to denote a number close 
but strictly bigger/smaller that 7. Having all these in mind let us gather the 
main results regarding (IFPDl) in one single statement. 

Theorem 2.2. Let 6 > 0 be fixed. Let (uo,Wo) be a pair of nonnegative func¬ 
tions in L^ n (H). Then, there exists a unique nonnegative weak solution to 
(IFPDl) . In the case i5 = 0 uniqueness continues holding up to a constant in the 
pheromone p distribution. Furthermore, the following estimates are satisfied by 
the solutions for any 0 < t < T <00; 
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(i) L'^-integrability 

IkWIU + II^WIU < C{mo,j)(l+ ^(i/!y/)+ ) : for any 7 G [l,oo), 

(ii) L°°-bound 

lk(i)||oo + \\u{t)\\^ < C{mo)(l + . 

The constants in (i) and (ii) depend on (El) but do not depend on T. 

(Hi) If, in addition, the initial data are in for 7 G (l,oo], then the 

previous estimates are improved to 

||w(t)||^+ + ||u(i)||^ < C'{mo,-f), 

where the constant C depends on El and on ||wo||^+, IIM 0 II 7 but not on T. 

Remark 2.3. Once the L°“ estimate has been established, it can be used to 
investigate further regularity of the solution. In particular, it implies that Vu 
and Vp are bounded on x O', for any 0 < t* < T and any domain O' 

strictly included in O; see [251 Th. VII.6.1]. Assuming that N,c,v are (7°°, we 
can boil down a bootstrap argument to show that the solution (m, w ) is actually 
of class C°° in any such subdomain [t^,,T] x O', see for instance [151 Prop. A.l 
& Th. A. 1]. 


3 Analysis of the model flFPD 


In this section we provide a series of a priori estimates to build up enough 
regularity to prove the existence of classical solutions of (IFPDI) . 

Definition 3.1. A classical solution {u,w,p) of the system (IFPDI) is defined 
as a triple {u,w,p) satisfying the following: 

(i) The triple (u, w,p) G (^([O, Tj; L^(0)) and each of the terms in the system 
(iFPDl) (i.e. dtu, Am, V • {uVp), and so forth) are well defined functions 
in L2 ((o,T) X 0), 

(ii) The equations (IFPDI) are satisfied almost everywhere, 

(Hi) The initial data = {uo,Wo) and the boundary condition (ED 

are satisfied almost everywhere. 


3.1 Prom to U* regularity. 

In this section we prove instantaneous generation of L'^ integrability, with 7 > 1 , 
for nonnegative classical solutions of (iFPDl) associated to initial data lying in 
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only. A crucial fact used in the argument is mass conservation, that is, 
classical solutions (u{t),w{t)) satisfy (11.41) that we rewrite as 


/ u{t,x)dx+ / w{t,x)dx= / Uo{x)dx+ / Wo{x)dx = mo for any t > 0. 

t/Q ^ Cl 


Proposition 3.2. Let {u,w) be a classical nonnegative solution of the system 
(IFPDI) with boundary conditions dni). Then, for any 7 G [l,oo) we have the 
estimate 

IkWIU + IkWIU < C'("lo,7)(l + ^ > 0, 

where the constant C{mo, j) depends additionally on (iHl) . but, is independent of 
time. Furthermore, for any 7 G [l,oo) this estimate can be upgraded, provided 
we add the dependence of the integrability of the initial data to the constant, to 


I (t) dx + f u'^{t)dx< C {mo, \\uo\\-y,\\wo\\j+), t>0. 

Jn Jn 

The constant depends on (iHl) (but is independent of time). 

Proof. We start by multiplying equation (iFPDil by It'*' with 7 > 1 , and 
integrating with respect to the space variable. We obtain 


f u'^dx-'^—- 

7 own 7 


f Vp • V('u'’') dx + (7 — 1) f ^|Vupdx 

Jq Jn 


= — u'^cdxF / Nwu'^ ^ dx . 


Now multiply the third (pheromone) equation by and integrate by parts to 
conclude that 


Next using 


/ Vp-V(u^)dx< / wu'* dx . 

Jn Jn 

[ dx = [ |Vu'’'/^pdx, 

Jn 7 Jn 


and the Young inequality with conjugate exponents p = 'y, p' = and 

p = 7 + 1 , p' = , respectively, we find 


— / u'*' dx + 4 
dt 


In 


7-1 

7 Jn 


|Vu'>'/2pdx 


< 7 / Nwu'^ ^ dx + (7 — 1) / wu'^ dx 
Jn Jn 

^lloo [ dx + 7 ||A^||oo [ u'^ dx 
Jn Jn 

+ i f dx + 76 ^/'’' f u'’'+^dx. 

£ Jn Jn 


(3.1) 









We have also used weighted Young’s inequality 


ab = 


.1//3 




/S' 


with free parameter e > 0, in the last inequality. A similar procedure on the 
second equation in (iFPDl) gives, for any a > 1, 


— / re" dx + 4 
dt 


a — 1 


In 


Jn 


|Vw“/2pdx 


<a cuw°‘ ^ dx + a(a — 1) / w°‘ ^Vw-Vudx. 

Jq Jq 

The first integral of the right hand side is directly estimated by 


dx + a / u;“ dx I , 


llcllc 


as a consequence of Holder and Young inequalities. There are two ways to 
estimate the last integral, depending on the assumptions on v: 

• We recognize aw°‘~^Vw = x 2Vw“/^ and we simply use the Cauchy- 
Schwarz inequality to obtain 


w°‘ ^Vw-Vudx 


Y{a- 1) [ 

Jn 

< ||Vu||oo ^ J |Vw“/^pdx^ ^a(Q! — 1) J dx 


N 1/2 

V 


which leads to 


— / w°‘dx + 2 


a — 1 


dt 


In 

< Hell 


o: Jn 


|Vu;“/2pdx 


(3.2) 


^ u“ dx + a(||c||oo + ||| Vu||^) w°‘ dx. 

When (d holds, we can integrate by parts so that 

a[a—l) f w°‘~^Vw-Vvdx = {a — 1) j Vu>“ • Vu < (a — 1) f w°‘Avdx 
J Q, J Q, Jci 


which yields 

^ f 

dt 


re" dx + 4 


a — 1 


'n 
< llcl 


|Vw“/2pdx 


||c||oo / M“dx + a(||c||oo + ||Au||oo) [ w“dx. 
Jn Jn 


(3.3) 
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Estimate (I3.2|) is better in terms of required regularity for v, however, (13.31) will 
play a role when a —>■ oo. The latter will be used to show uniform propagation 
of the L°°-norm later on. 

Now, we use the fact that the space dimension is n = 2 and we appeal to the 
following Gagliardo-Nirenberg-Sobolev inequality (see e. g. [Jl eq. (85) p. 195]), 
which holds for any a > 1 , 


CQ!+1 


dx < 


< 


c(o,a)iieiii 

C{n,a) f ^dx( f rd 2 ;+ f |V(r^2)|^da; 

J Q, \J Q Jq 


(3.4) 


Combined with J u(t, x) dx+f iu(t, x) dx = mo, it allows us to absorb the higher 
exponent of u in the right side of by choosing e > 0 sufficiently small. We 
can thus find positive constants C(mo) and C'{mo), depending on (iHl) . and on 
the exponents a and 7 , such that 


_d 

dt 


f u'^dx + C'{mo) [ dx 

JQ Jn 


< C{mo) 


V? dx - 


/ re'*' dx + 

Jn 


J dx^ , 


(3.5) 


_d 

dt 


w°‘ dx + C{mo) 



dx < C{mo) 


dx + 



In order to control the right hand side of inequalities (13.51) exponents in the left 
hand side have to be bigger than exponents in the right; therefore, we are forced 
to choose 7 < q; < 7 + 1. We shall use Lebesgue’s interpolation inequalities 



01 - e ( 0 , 1 ), 


»2 = 5gSie(0,i). 

ikiu<Nir'iK^+i, 

03- e(0,i). 

Let US introduce 


U{t) '■= [ yJ{t)dx, 

W{t) := [ w°‘{t)dx. 

Jn 

Jn 

>From (|3.5I) we are thus led to 


±U + C'{mo)\\u\^X\<C{mo) 


Aw + C'{mo)\\w\Cll<C{mo) 

(vE + ihii5;.Y^®^). 
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For any e > 0,/3 > 1, we can find C{e,P) > 0 such that s < es^ + C{e,f3) for 
any s > 0. Thus, adding up these inequalities and setting Z := U + W, it follows 
that there exists constants Ci{Tno) and ( 72 (mo) such that 

^Z + Ci(mo)(||w||:(+} + |lHI“i}) <C 2 (mo)(Z + l). (3.6) 

Using Lebesgue’s interpolation again 

we readily infer 

+ ikiiail > c(mo)(hi|("+i)/^^ + 

> C'(mo) > C{mo) — C{mo) , 

by using the fact that s i—>■ is non increasing. Using s < es^ + C{e, j3) again, 

and coming back to (13.61) . we arrive at 

4 ^ + C((mo)Z^ <C'(mo), 

df 

with constants C{{mo) and C^imo) depending on the initial mass mo, ((Hi) and 
the exponents a, 7 . Therefore, the comparison principle in Corollary I A. 2 1 yields 

^(t)<C(mo)(l + ^). 

In other words, for any 1<7<Q!<7+1< 00 , we have 
||w(t)||a < C{mo)(l + , 

||M(i)|l7 < C{mo)(l + -4r) = Cimo)(l + ^4)+) ' 

In the last equality we have used the fact that 7 < a can be taken as close as 
desired. Furthermore, by invoking Lemma lA.ll the estimate on t 1 —> Z{t) can 
be upgraded to supoQ Z{t) < C provided we add the dependence on the norms 
||uio||a and IIU 0 II 7 in the constant. This completes the proof of Proposition 13.21 

□ 


04 = 


05 = 


72-1 
7^ ’ 

a2_l 


a 


2 ’ 
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3.2 Prom to L°° regularity. 

In this section we prove further gain of boundedness for (IFPDI) . showing that 
classical solutions are in fact bounded for any positive time. We adopt De 
Giorgi’s energy method m, which has been successfully adapted to investigate 
the regularity of PDEs systems see e. g. EllZllli]. Let us start the discussion 
and we shall state the theorem at the end of this section. 


Consider a classical non negative solution w of the equation 
dtw — Aw + V • {w\7v) = f , in [0, T] x 17 , 


(3.7) 


with V and / given functions. For the boundary condition we assume that 
(Vru — wWv) • n = 0. Define the level sets 

WA = (w - A) , A>0. 

Multiply (13.7|) by w\. Owing to the boundary conditions, we obtain 

— / dx + 2 ( |Vwa|^ dx <2 f wVv ■ Vwx dx + 2 f f+w\ dec. 

“ Jn Jn Jn Jn 

Young’s inequality leads to 

■ / w^dx+ |VuiA|^dx< / |uiVn|^ 1{^>A} dec + 2 / f^w\ dx. (3.8) 

jn Jq Jci j o. 


_d 

dt 


We split the reasoning into two steps. We start by proving the L°° bound on 
a given time interval [t*, T\. Secondly, we shall extend the estimate to infinitely 
large time intervals. Thus, let us consider 0 < < T < oo. Let M > 0 and 

define the following sequence of levels and times 

Afc= (1-1/2'=)M, 4 = (l-l/2'=+i)G, fc = 0,1,2,... 

Define the following energy functional for the level sets 

Wk ■= sup f wl dx + ( f |Vwfc|^da;df, (3.9) 

te[tk,T]jQ Jtk Jn 

where we adopted the notation Wk '■= w\^. With A = Afc, let us integrate 
inequality (13.81) over the time interval [s,t]; we obtain 

/ w^{t,x)dx+ / / IVwfe 1^ da;dt'< / w^{s,x)dx 

Jq J s J q J q, 

+ f f |wV-c|^ dxdt'+ 2 f f f+Wkdxdt'. 

s ^ ^ s ^ ^*2 
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We use this relation with tk-i < s < tk < t < T. It implies 


Wk < / wf{s,x)dx 

Jn 

+ / |wVu|^ dxdt'+ 2 / / f+Wkdxdt'. 


We take the mean over s G t^], bearing in mind that tk — tk-i = 

It yields 


2/c+l pT p 

Wk < - / w^dx ds 

J tk_i J n 


+ 



wVv 


d{u;>Afc}d^dt +2 



/ /+ Wk dx dt'. 

Jn 


(3.10) 


We are going to make use of the Gagliardo-Nirenberg interpolation inequality, 
see O eq. (85), p. 195], 

||w||]( < C'(n,p,a)||w||^“ 1 = Q - (3.11) 

which holds for any a G [0,1] and 1 < p, </ < oo (note that we perform the esti¬ 
mates without restricting the space dimension for the moment). Thus, choosing 
ap = 2 it follows that 

||u;|]^ < C(II,n)|]u;f^.|]u;||r^ , p = 2^ . (3.12) 

Note that if rcfe > 0 then Wk-i > 2~^M, and, as a consequence, 

1{»>A,} < (y , Va>0. (3.13) 

Having this in mind one can play with the homogeneity of the right hand terms 
in the level set energy inequality (13.101) . Indeed, with a = 4/n, the first of them 
can be evaluated as follows 


2^+1 pT 


f 2 ^^^ f / 2 ^ \ ^ 

y^u;^l{„>A,}da:ds < y— ^ j^wl_^i^—Wk-iy dxds 


/ i fe _ 1n 

/.T 

< 2 


Jtk-i Jn 


Wk-\ dx ds 


< 2C'(H,n) 


l-T 


M^t 


★ ^tk-1 


||wfe-l||2 + ||V?Ufe-l||2)||wfe-l|]2 " ^ds 


2^^= 21±2 

< 2C(H,n)(l + T)^^W, 




k—1 ’ 


(3.14) 
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by virtue of (13.121) and the definition of Wk-i- The last two terms in the right 
hand side of (13.101) can be treated by using a similar procedure, together with 
the application of Holder’s inequality. On the one hand, bearing in mind that 
the tfe’s are all larger than t*/2, we get using p.l3l) with a = p 


f I 

Itk-i JQ 


wVu|^l{„>Aj,}dxdt < J \\wVv\\lg,(^J ’ dt 

pT 




< sup \\wVv\\lg, 
Ml t>i^ 


— 1 


n 




< C{^,p,a)i 


h 

~ 


sup ||wVu|| 2 ,- 





tk-1 


We have used (j3.11|) in the last inequality. We choose the parameters so that 
= 2 for some 0 < a < 1, which can be achieved as long as 1 < q < (more 
precisely, going back to (|3.11|) . note that a = ij^^qjn p = 2(1 + 2g/n)). It 
follows that 


/ / |■u;V^;|^l{u,>A^,}da::d^ 

J tk — 1 ^ 

< {l+T)C{il,p,a)^^ (supt>^ \\wVv\\lg,'j . 

On the other hand, since Wk < w have 


(3.15) 



/ /+ Wk dx dt < 
Jn 


[ sup ||w/+|| 
\‘>¥ 








dt 


^ 77 e ( sup ||w/+||q. ] f ( f wl_^ dx) ’ dt (3.16) 

i2^^ [ \ 1. 

< (1+ T)C7(0,p,a)i^ sup ||/+|| 2 ,'|k|| 29 ' ITfc“_i, 

V‘>¥ / 


still using the same relation between p, g and a. We go back to (13.101) . with 
(I3.14|) . (13.151) and (|3.16l) : we arrive at 


Wk <{l+T) 2C{n,n) 


4+n li 

^ + 2 


M^t. 


-Wk-, 


+ 20(0, p, a) 9—^ sup (||Vu|l 


U\M\lq' + \\f+hq'M\2q)wi^ 


(3.17) 


For the final step, we specialize to the case of space dimension n = 2. As a 
matter of fact, we observe that 1 < g < = oo. We can then appeal to the 
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estimates in Proposition 13.21 which imply (particularizing to /+ = m c) 


sup (||V?;|l^|lu;||^^, 

t>-^ 


||/+||2g'|kl|2g') < C(TOo) (l + • 


Therefore, (13.171) becomes 


Wk < C{1 + T) 



wLi 


+ 




(^)+ 


2(9+1) , 


M 


2(9+1) 



(3.18) 


(3.19) 


We can take advantage of the fact that the powers = 2 and are strictly 
larger than 1 to conclude. Indeed, a direct calculation shows that Wq is 
a super solution of the first order difference equation (13.191) with choices of 
a G (0,1) small, and then M large (see e. g. |151 Lemma 3.3 & 3.4] for similar 
arguments) such that 

2(9+1) I'l 

max < 2'^a, 2 i a® > < 1, 

r CWo CWn’ 1 M 1 

I a2 MH^ ’ ^2s±i (s^)+J| - 2(1+ T)' 

6* 


By a comparison principle, we check that Wq > Wk holds under these cir¬ 
cumstances. We conclude that 


lim Wk = 0 . 

k—¥oc 


Let us suppose 0 < t* <C 1. Then, observe that such a choice of M explicitly 
takes the form 


M = max ■ 


l2C{l+T)Wo J {2C{l + T)) — Wo^ 


aHi, 


( 1 / 2 )- 


t 


-)}, (3.20) 


where the constant C depends on mo, and of the exponent q. We clearly 
have ^ 

Wk> / / w^(t,a:) l{,„(t,9,)>M(i-i/2‘^)}da;df. 

J t* Q 


/ ti, JQ 

Letting k run to 0, by virtue of Fatou’s lemma we deduce that 

-T 

T-. 


—r I I '^{w{t,x)>M}dxdt = 0, 

— t* Jt„ Jn 


which eventually implies 

0 < w(t,x) < M a. e. x 17. 

This proves most of the following result. 


(3.21) 
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Proposition 3.3. Let {u, w) be a classical nonnegative solution of (IFPD|) with 
boundary conditions (CH). Then, the following L°°-estimate holds 

||u(t)||oo+|kW||oo<C(mo)(l + ^), tG(0,T], (3.22) 

where the constant C depends additionally on IHI, but, it is independent of 
T > 0. Furthermore, when (HU holds, estimate (13.221) can be upgraded by 
adding the dependence on the -norms of the initial data in the constant, 

max{||'u;(t)||oo, ||w(<)lloo} < ||wo||oo, ||mo||oo) , t>0, (3.23) 


where the constant is still independent of the time T > 0. 


Proof. Proceeding as we did in (j3.8p . we can establish the following inequality 
for w: 


i 


da; + / \Vw\ dx< kcVf da: + 2 / Z+wda:. 


(3.24) 


We integrate over [s, f\ and we get, with /_|_ = uc, 

/ w^{t,x)dx-\- / / |V?c((j,a:)|^ da:dCT 

J Q. J s J Q 

< / w^(s, a:) da: + II Vw||oo / / w^{a,x)dxda 
J Q, J s J O, 

+ l|c||oo / / + w^)(cr, a:) da:dcr. 

Js Jn 

We use this relation with s = t*/2 < t < T, which allows us to obtain the 
following estimate (recall the dehiiition of Wq in (|3.9I) 1 


W'o < \\w{U/2)\\ 


'(l|Vu||oo + ||c||oo)(r-^)( sup ||u:(t)||2+ sup ||u(t)||^) . 


By Proposition 13.21 Wq is thus dominated by C'(mo)(l + T)(l + 1/tf^). 


(3.25) 


Now, we set T = 1, say. Going back to (|3.20l) and (13.211) . we conclude that 


||w(t)||oo < M < C{mo) -Y+ 


(3.26) 


holds for short times, say 0 < t* < t < T = 1. It is clear that the same reasoning 
can be applied for any T, and it would provide a bound depending on T. In order 
to justify that the estimate can be made uniform with respect to T, we thus 
proceed differently, by extending the estimate obtained on [t*,l]. Indeed, for 
any ti > 0, the shifted function (t, x) i—>■ {t, x) = w(t-l-ti, x) is still a solution 
of equation (|3.7p . with data wt^ (0, x) = w(ti , x) and the appropriate right hand 
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side /. In particular mass conservation implies f^(u(ti,x) + w(ti,x)) dx = mg, 
so that the constant C(mo) does not change. Therefore, we pick 0 < ti < T 
and we can repeat the same arguments for , which leads to the same 
estimate (I3.26|) as for w on the time interval that is to say (13.261) holds 

for t G [ti + + T], We have thus extended (13.261) over + T], see 

Figure m and we can repeat this procedure. This completes the proof of p.22p 
for w. 

IkWIloo < 


0 u 


ti ti +T ti -\-T 

ll'^(i)l|oo < C{mo)/t\* 

Figure 1: Extension of the local result. 


Next, using elliptic regularity for the pheromone equation, see [T], and 
Sobolev’s embedding (still with space dimension n = 2), it follows that 

||Vp(t)||oo < C\\p{t)\\^2,2+ < C{n) ||w(t)||2+ < C{mo) (l + ■ 

Equipped with this estimate we repeat for u the arguments used for esti¬ 
mating w, just changing c to N and replacing the uniform estimate on Vu 
by this time-dependent estimate for Vp. This finishes the proof of estimate 
(I3.22p . As a matter of fact, note that (13.261) applies in particular to t = t*: 
||w'(^*)||oo < C{mo)/tl^ holds for any 0 < t* < 1. 


For proving the uniform bound (13.231) let us set 


A := max 


I sup||u(s)|| 

s>0 


2 +) 


sup ||iy(s) 

s>0 


I2+ 1 < 00. 


(3.27) 


We slightly modify the previous analysis. In particular, we change the definition 
of the tfc’s by setting tk = t*(l — 1/2^), so that it now starts from tg = 0. Going 
back to (13.241) . we get Wg < Urcolli + CT, with C depending on (IHl) and the 
constant A in (13.271) . Furthermore, the additional information in (13.271) allows 
us to control uniformly with respect to time the norms (bearing in mind 
that making q large means q' close to 1) so that (I3.18P becomes 

sup (||Vx||^||u;||ig, -k \\f+\\2q'\\w\\2q'^ < C. 


Accordingly, inequality (13.191) is changed into 


Wk<Cg 


23k 




,As±llk 


-W, ", 

2(g+l) k — 1 ’ 


M- 
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where Cq and Ci depend on but also on T and A, while they do not 

depend on k and t*. Reproducing the same argumentation as above permits us 
to establish that 

||w(t)||oo < , (3.28) 

holds for 0 < t < 1. Furthermore, elliptic regularity [T] yields 


||Vp(t)||oo < C(n)||«;(t)|| 2 + < A (3.29) 

and, as a consequence, we get a similar estimate for u: 

\\uit)\\oo<C{mo,A)(l + ^'^ (3.30) 

holds for 0 < t < 1, too. Let us set 

U{t) := ( u'^{t,x)dx, W{t) := f w'^{t,x)dx. 

Jq Jq 

Since (O holds, we can make use of inequalities (EH) and (13.3p (with 0 = 7 ), 
and we get 

^ < ||A^IUt^ + 7(l|iV||oo + lk(0lloo)i7 

< ||A^||ooLF + 7 r(i + ^)i7, (3.31) 

^<I|c||ooC/ + 7(I|c||oo + ||ahU)ir, 

for t G (0,1), with B > 0 independent of 7 > 0. Adding the two inequalities in 
(|3.3ip . we deduce that Z := U + W satisfies 


dZ 

- < 

dt 



Z, 


for t G ( 0 , 1 ), and a certain constant C > 0 which does not depend on 7 . 
leads to 


Z{t) < Zoexp (^ 7 ^ (1 + ds) < 


which eventually implies the estimate 


It 


max{||'u(t)|l.^, 11 ^( 011 ^.} < max{||Mo||.^, 


which thus holds a.e. t G (0,1). Sending 7 —?> 00 proves estimate (I3.23p for short 
times. 

Going back to the second part of the statement in Proposition [321 we realize 
that (13.271) is ensured when the initial data belong to (7 This obser¬ 

vation ends the proof of the L°“ estimate for 0 < t < 1. Combined to ([3.2211 . it 
provides a uniform estimate on the solution for any t > 0 . □ 
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3.3 Construction of classical solutions 


We have now the tools to prove the global well-posedness of the system (IFPDI) . 
Firstly, we establish the existence-uniqueness of solutions locally in time, and, 
secondly, we extend the result, as a consequence of the obtained estimates. We 
start by working with smooth and bounded initial data, say (uo,Wo) S C^{Q). 

Theorem 3.4 (Classical solutions). Let Uo-, Wo G be nonnegative initial 

data. Then, for every T > 0 the system (iFPDl) with boundary conditions (HI) 
admits a unique nonnegative classical solution {u,w,p) S y which satisfies the 
estimates of Propositions \3.2\ and \3.tA 

Proof As explained above the proof splits into two steps. 

Step 1: Local existence. Let us introduce the convex set 

y = {e e C{[0,T]■L\VL))c^L^{0,T■H\n)), 0 < f{t,x) < 2||u;„||oo}. 

This space is endowed with the norm 

Ikb = sup \\u{t,-)\\L 2 . 

0<t<T 

Let ip G y he such that (p(0) = Wq. To such ip we associate the solution p{p) of 
the pheromone equation 

—Ap = p — Sp, (3.32) 

and, then, we associate the solution u{p) = u{p,p{p)) of the equation 

dtu — Au + V-(^uyp{p)) =—UC + pN, u{0) = Uo . (3.33) 

Finally, let w{p) = w{p,u{p)) be the solution of 

dtw — Aw+ \7-[wVv) = uc —wN, w{0) = Wo . (3.34) 

Equations are complemented with the zero-flux boundary conditions (11.11) . Note 
that (|3.32|) . (|3.33l) and (13.341) are linear equations and existence of solutions is 
provided by the standard theory, see e. g. 0 Theorem X.9]. In particular both 
u{p) and w{p) belong to C([0, T]; L^(n)) fl L^(0, T; iL^(n)). We aim at finding 
T small enough, so that solutions of (IfTO can be obtained as fixed points of 
the mapping 

$ : p I —> w{p). 

We shall prove that d* is a contraction in y, for a sufficiently small time T > 0 
depending only on m and the initial data (uo,Wo)- This is a consequence of 
the maximum principle for the linear equation 

dtif — Atjj + V • (Btp) +b'tjj = f, 

complemented with the boundary condition 

V^/> • n — fiB • n = 0. 


19 














We assume (related to (11.21) 1 


^>0, />0, S-n<0. 


Then, the solution ip associated to ip{0,x) = ipo{x) > 0 is non negative. This 
result can be obtained by considering the identity 


/ G{ip)dx-\- [ G" {ip)\V ipl"^ dx 

dt Jn Jn 

= [ ipG"{ip)B -Vipdx- [ biPG'{ip)dx+ f fG'{ip)dx. 
J ^ J ^ t/ 


\z]^ 

We work with the convex function G{z) = [z]- = min(z,0) < 0. It follows 

that G'(z) = [z]- < 0, and zG'{z) > 0, zG"{z) = G'{z) so that \ipG"{ip)Vip\ = 
y/2G{ip) X y/\G"{ip)\Vip\. Then Cauchy-Schwarz and Young inequality yield 


_d 

dt 


[ G{iP)dx<\\B\\l [ G{iP)dx. 

J Q. JQ 


We conclude as a consequence of the Gronwall lemma that ip{t,x) > 0. Let us 
set 

Observe that limt_>o7(t) = ||'0o||oo- Applying the previous maximum principle 
to 'tp{t,x) := 7 (t) — ip{t,x), we conclude that ip{t,x) < ^{t). Indeed, note that 
ip(t,x) satisfies the non homogeneous Robin’s condition {y-ip — Bip)-'n. = -jB-n 
with a non negative source —^B • n > 0 so that the computation now involves 
the boundary term — B n jG' (j — ip) da which contributes negatively owing 
to the orientation of B towards the interior of fl. 

Since 0 < <p{t,x) < 2||wo||oo holds on [0,T] x O, we have that p and Ap are 
bounded in L°°{{0,T) x fl), and Vp as well. This allows us to use the above 
maximum principle with the solution of (j3.33|l . and then also to the solution 
of (j3.34p . to justify that 0 < w{t,x) < 2||wo||oo holds a. e. {t,x) G [0,T] x fl, 
provided T is small enough. We thus have d'Cy) C y. 

Let us show that d) is a contraction in y, possibly at the price of making 
T > 0 smaller. Let pi, p 2 & 'P and let (ui, Wi, pi), with i = 1,2, be the 
corresponding solutions of respectively, p.32|) . (I3.33|) . and p.34l) . Define w := 
wi — W 2 = 4>((/3i) — 4)(<P2) and so forth. We compute 


l_d 
2 dt 


wP dx 


iVwf dx - 


N\ 


wP dx 


= / cuw dx 

Jn 


w'S/v ■ Vw dx 




llcIlL) / 

Jn 


tep dx + i 
n 2 


dx - 


luP dx. 


Gronwall’s lemma thus yields 


\w\^ 


dx < e 


Kit 


luP dx ds 


0 Jn 


(3.35) 
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with Ki = ||Vw||^ + ||c|l^. We proceed similarly with (13.331) and we get 



where we have used the equation dpi — Api = wi. Since 0 < x) < 2||r(;o||oo 
holds on [0,T] x fl, we deduce that p 2 and Ap 2 are bounded in L°°((0,T) x fl), 
and thus, the maximum principle applies for (j3.33ll and we can assume that 
||rt 2 ||oo < 2||uo||oo holds on [0,T]. Furthermore, we have ||Vp ||2 < with 

C depending on i5 > 0. It follows that 



holds on [0,T]. Hence, we get 



with K 2 = ||.^V|1^ + 2||ryo||oo. Finally (going back to (j3.35|) l. we find K,M > 0, 
depending on m and the L°° norms of the data (uq, Wq), such that 



holds on [0,T]. This justifies that d* is a contraction on y, for a small enough 
time interval. The fixed point defines a solution of (IFPDI) . 

Step 2: Global solutions. We remind the reader that we are working with data 
{uojWo) S (7“(n). The obtained solution is bounded on (0,T) x H. Therefore, 
the classical theory of parabolic equations tells us that Vu and Vp are bounded 
on [t^,, T] X n', for any 0 < t* < T and fl' strictly included in 12; see for instance 
PSI Th. VII.6.1]. We can then perform a bootstrap argument to show that 
the solution is C°“([t*,T] x H'), if the coefficients N,c,v are assumed to be 
C°°. Moreover, by virtue of {uo,Wo) G C^{i2), the data satisfy compatibility 
conditions with the boundary conditions, so that we have actually constructed 
classical solutions, see [261 Th. 4.3] or [25l Th. V.7.4]. In particular, the 
a priori estimates discussed above apply. Proposition 13.31 provides a uniform 
L°° estimate on the solution, which depends only on © and |]uo|]oo, |lwo|]oo- 
Therefore, we can reproduce the fixed point argument and extend the solution 
over [0,oo). □ 
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We finish this section studying the stability of classical solutions in the space 

§ := I {u, w) G (0, T; n (n)) } , T > 0 . 

Proposition 3.5. Nonnegative classical solutions of (IFPDI) with boundary 
conditions (EH) are stable in S. More precisely, given two initial data (ui^o, Wi,o) G 
with i = 1,2, then solutions {ui,Wi) satisfy the estimate 

||Ml(t) - U2{t)\\2 + 11^1 (t) - W2{t)\\2 

< (^||mi,o - U2,o\\2 + Ikl.o - , t>0. 

The constant C{mo) depends on (iHl) and the || • II 2 +-norm of the data. 

Proof. Let u = ui — U 2 and similarly for the other differences. The system for 
the differences reads 

diu — Au + V • (uVpi) + V • {u 2 S/p) = —cu + Nw, 

dtW — Aw + V • (uJVr) = cu — Nw , (3.36) 


—Ap = w — 6p. 

Multiplying the second equation in p.36l) by w, integrating by parts and using 
Cauchy-Schwarz, we are led to 

~\\w\\l+\\yw\\l < |N|2(||Vz;||oo||VuJ||2 + IICII00INI2) , 

which yields 

^|H|2 + ||Vw||^<c(|H|2 + ||u||2). (3.37) 

for C = llVull^ -I- ||c||^. Next, using the equation for u we get 

+ llVull^ < 2 (||uVpi||^ + \\u2ml) + INI2 + , 

The right side terms above are easily controlled using (13.2711 . (j3.28|) . (13.2911 . 

(iTiinll 

\\umpiit)r < \\ypimi\\m\\l < cwmwi 

||u2(t)Vp(t)||^ < \\u2m%,\\ym\\l < (ll«2(t)||}/^V2(t)llL+'/^)l|p(i)ll?^2 

- ^(1 +<7e[l,oo). 
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The constant C now depends on nrio and on the norm of the data. Taking 
q arbitrarily large, we obtain 


+ ^l|V^ll2 < c[l + {\\u\\l + INI 2 ) • (3.38) 

Now define Z{t) := ||m||| + ||uJ|| 2 . Thus, adding up estimates (13.371) and (|3.38|) 
it follows that 


dZjt) 

dt 




Integration of this ODE leads to the result. 


□ 


3.4 Global well-posedness of weak solutions 

We are going to prove the global existence of weak solutions, in the sense of 
Definition 12.11 

Theorem 3.6 (Global well-posedness). Let arbitrary T > 0, 6 > 0 be fixed and 
assume nonnegative initial data {uo,Wo) & nL^^(D). Then, there exists a 

unique nonnegative weak solution for the system (iFPDll . Such solution satisfies 
the estimates of Provositions \S.2\ and lS.tA 

In the case ^ = 0 uniqueness continues holding up to a constant in the 
pheromone p distribution. 

Proof. Take a sequence of non negative initial data (u*, w^) G (^“(D) such that 

(mq, w^) (uo, Wo) strongly in O (D) . 

Using Theorem 13.41 on global well-posedness of classical solutions, we have a 
sequence [u^,p^) G C(0,T; L^{Ll)) of solutions to the system (jFPDI) . It 
is not difficult to check that, in fact, such a sequence is uniformly bounded in 
L'^ (0,T; H^{fl)) with time derivatives uniformly bounded in L^(0,T;i7^(D)*). 
Let 

Zkiit) = Wu’^it) - u\t)\\l + ||n;'=(t) - wHt)\\l , k, I > 1, 
be the Cauchy differences. Note that Proposition 13.21 implies that G 

L°° (0,T; (D)) uniformly with respect to k > 1, hence, the stability result 

of Proposition 13.51 implies that 

Zki{t) < Zki{0)e ^*, 

with constant C independent of the indices k, 1. In this way we conclude that 
both and are Cauchy in L°° (O, T ; L^(D)). Thus, the following convergence 
properties hold: 

w^) {u, w) strongly in L°° (O, T; ^^(D)) , 

(u'^,w^) ^ (^u,w) weakly in L^(0,T;i7^(D)) 

{dtu^,dtw^) —>■ {dtu,dtw) weakly in (O, T; i7^(D)*) . 
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With standard estimates for the elliptic problems, we deduce that is bounded 
in L°°(0, T; and Vpfe converges to Vp strongly in L°“(0, T; L^(n)). Ac¬ 

cordingly the product UkVpk converges to uVp strongly in b^((0, T) x ft). As a 
consequence, 0 < {u,w) G L°°(0, T; L^(n)) is a weak solution of (IFPDI) . Note 
that the initial condition (u(0),rc(0)) = (uo,Wo) is satisfied by continuity at 
t = 0 which follows from the estimate on the time derivatives. 

Finally, using approximation by classical solutions again one proves that the 
stability result of Proposition [3^ is valid for weak solutions as well. Uniqueness 
follows from here. □ 

Remark 3.7. It is likely that the integrability of the data is not optimal 
for the theory of existence of solutions. However, it simplifies the analysis in 
two directions. On the one hand, it makes the dehnition, and the stability, of 
the product uVp meaningful. On the other hand, the underlying estimate is 
also useful in the proof of Proposition 13.51 which implies the uniqueness of weak 
solutions. 


4 Analysis of the model fISPDI) 


In this section we perform the analysis for the system (ISPDt that includes 
regularity estimates and global well-posedness in the case where the problem is 
set on the entire space O = K.^. A short comment for the case O = R.” with 
u > 3 is included in the Remark 14.61 Details are given only in the arguments 
for obtaining a priori estimates since the ideas to make such estimates rigorous 
were previously presented (and are quite classical). We remind the reader that 
the main hypotheses imposed on the parameters are gathered in m- Eq. 
(Ispnli involves an evolution equation for the food concentration c; therefore, 
in addition we will assume 


Co G L°°(K") , P G L°°(R"). (H+) 

Of course, we continue assuming that the total population of ants is integrable: 
Uo and Wo belong to and we still denote 

/ {uo + Wo){x) dx = mo- 


4.1 Estimates for solutions of the heat equation 

The strategy for proving estimates for the system (ISPDt is based on a sharp 
control of the gradient of the pheromone Vp in terms of the ant population w: 
recall that for the system (iFPDt such an estimate was a direct consequence of 
the regularity analysis for elliptic equations, a la Agmon-Douglis-Nirenberg [1]. 
We substitute this argument with a direct analysis of the Duhamel formula for 
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the heat equation. To be more specific, let (p be a solution of the problem 
dt^ — Aip = f in (0, oo) x R", 




(4.1) 


t=o 


= ipo on . 


Then, (p is given by the explicit formula 


ip{t,x)=f K{t,x-y)ipoiy)dy+ f f K{t - s,x - y)f{s,y)dyds (4.2) 
Jr" Jo Jr" 


where 


1 _l£li 

K{t,x) = ——e 


as long as it makes sense. Throughout this section the technical lemmas will be 
presented as general results valid for any space dimension n> 2. 

Lemma 4.1. Fix q G (n, oo] and 9 G (0, 2] suck that 

Oq 


^>1 

2 ’ 2-9 


> n . 


Additionally, assume that for some T > 0 

Vifo G L«(R"), f G L°°(0,T;LinL‘?®/^(K”)) . 

Let p he a solution of the heat equation (14.11) . Then, for any t G (0,T] 

||Vv5(t,-)||,< ||V(^o|l, 

+ Cn (1 "b 


^(n-q- 


In/ \ e q{n-\)-r, 

t _Vi I iir/„.Mi.i ..... iir/. .Ml" 

q')/2q’ )y 


1 + sup ||/(s,-)||i) sup ||/(s,-)||,"e^ 

Q<s<t ^ 


(4.3) 


The constant Cn depends on the dimension n, q and 9. Furthermore, the supre- 
mum in time can be estimated for any T > 1 as 


sup ||Vv3(t,-)|l9 < 2||V(po||q 

0<t<T 


g(n-l)-n 

dq-2 


+ Cn(l+ sup ||/(s,-)||l) sup ||/(s,-)||"e /2 

N 0<s<t ' 0 <b<T ^ ' 

The constant Cn depends on the dimension n, q and 9. 

Proof. Using the explicit solution of the heat equation (14.2L we obtain 

Vip{t,x)= \JxK{t,x - y)Lpo{y)dy 

+ f f VxK{t - s,x - y)f{s,y) dyds =: Ii{t,x) + l2{t,x). 

Jo Jr^ 


(4.4) 
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Estimate of the term Ji simply follows from integration by parts and Young’s 
inequality for convolutions 


/ '^xK{t,x-y)(po(,y)dy 

[ K{t,x - y)W(poiy)dy < ||i4r(t, •)||i|| Vv^olU = || Vf/JolU- 


(4.5) 


Let us focus on the term I 2 . Fix 0 < e < t and consider the decomposition 

l2{t,x)= f f ...dyds+ f f ...dyds=:l2(t,x)+l2{t,x). 

Jo jB" Jt-e jB" 

In what follows, we will repeatedly make use of Young’s inequalities for convo¬ 
lution [3 Th. IV.30] 


\\f*9\\q < ll/lipbllr, 1/p+l/r = 1 l/q. 
For it yields 

|x-b|2 

-r 

\/(*-s) 

(x — y) 

r. AU- '' ^ J 


(4.6) 


-2(4i^/o " /b" e /(s,2/)d2/ds 


< 


< 


< 


2(47r 

1 


1_ /■*■" 

t)”/^ Jo 


(i _ s)(n-Hl)/2 

t—e Y 


I e 

B" \/(t — s) 

X _ , 


fis,y) dy 


ds 


2(47r)"/2 Jo {t - s)("+i)/2 

a:||^ 2g' supo<,<J|/(s, QUi 

4 "/ 29 ' 7 rn /2 ^ _ g/ ^{n-q')/2q' 


q \ 

da:j ||/(s,-)||i ds 


(4.7) 

Note that the case 9 = c» is valid. Next, we get (using (14.6|) with p = q912 and 
r = a) 


1^2 ' 




< 


< 


< 


2(47r)"/2 {t - s)("+i)/2 


/B" 

(x — 7/1 


B" 


V{t-s) 


f{s,y) dy 


ds 


2(4^)V2 (t - s)("+i)/2 Vir 
, 2a' 


■ 4(t-s) 


da: 


1 /ct 


)II 9^/2 ds 


/pn.l2a'^nl2 (j! _ ^ t-e<s<t 


sup ||/(s,-)||,e/2e^"'-”^/""' 


(4.8) 
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where a' = > n. (Note that the case 0 — 2, that is a' = oo, is a valid 

choice.) Thus, gathering (14.51) . (14.71) and (|4.8I) we are led to 


I^V3(ir)|L = \\h{t,-) + hit,- 


^ Ilv 7 II I /-> ll/('Si Olll (a'-n)/2cr' \\ n/■ \|i \ 

< ||V^o||^ + C„( - -^-_g,y2g’ - + ^ 

(4.9) 

with Cn depending on n, q, 9. We shall choose e S (0, |] defined by 


£ = min 


{l,|} 1+ sup ||/(s,-)llge /2 


2 


to equalize the homogeneity. Consequently, we obtain the estimate 


||Vv^(t,-)||,< liv^oll, 

+ ^(l + ^(173^)(l + ^supjl/(s,-)l|i) ^sup ||/(s,. 


a'in-q') 

n(cT'-q') 

qS/2 


where the constant C > 0 depends only on the dimension n, on the exponent q 
and 9. Estimate (14.31) follows by noticing that 


a'{n — q') q9 1 n — q/{q — l) 9 q{n — 1) — n 

nia'— q') 2 — 9n q9/[2 — 9) — q/[q — \) n 9q — 2 


Now, for estimate (14.41) fix T > 1 and gather (14.51) . (|4.7I) and (|4.8I) to obtain 


sup ||V(/?(t, •)|| = sup ||/i(t,-)+/ 2 (t,-) 

e<t<T ^ e<t<T 


— L “I" 


/ SUPo<s<T ll/(gr)l|] 

1 ^{n-q')/2q' 


.£(-'-n)/ 2 .' 3 ^p ||/(s ,.)||,./2 

0<s<T 


with Cn depending on n,q,9. Note that e G (0,T]. Similarly, when t G [0,£] it 
follows that 


sup ||V(^(t ,-)||9 = sup Whit,-) + hit, ■)\\q 

0 <t<£ 0 <t<£ 

(4.10) 

<W'^'-Po\\q + Cne^‘^' sup Wfis,-)\\qe/2- 

0<s<T 


Estimates f|4.9p and (|4.10p lead to 


sup ||V(/5(s,-)|l9 

0<s<T 


<2WWq^oWq + C, 


( suPo<^<tII/(s. 

I ^(n-q')/2q' 


-+e("'-")/2^' sup ||/(s, 

0<s<T 


\qe/2 


(4.11) 
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Again, we shall choose e G (0,1] C (0,T] defined by 


e= 1+ sup \\fis,-)\\ge /2 

\ 0<s<T , 


to equalize the homogeneity. Consequently, we obtain the estimate 
SUPo<s<T l|Vv5{s, Oil? 


(t' {n-q') 
n(<T'-g') 


< 2||V</5o||? + c[l + suPq< s<T I|/(s,-)||l) SUPo< S<T ii/(^r)ii;';;- 
where the constant C > 0 depends only on the dimension n, on the exponent q 


and 9. This proves (14.41) . 

Remark 4.2. The same result applies to the damped equation 
dtip — Aip = f — 5(p in (0, oo) x R", 




t=o 


= on 


□ 


(4.12) 


Indeed, it suffices to apply Lemma [4.II to e 

Lemma 4.3. Fix 

ipoGL°°{W^), Vu e L°°((0,T) X R") , 

/ e L“(0,T;Li(R")) (~lL°°((0,r) x R”) . 

Let ip be a nonnegative solution of the problem 

dtp> — Aip + V • {(fiVv) = / m R" X (0, oo) 

ip = ifo on R" X {t = 0} . 

Then, we have for any T > 1 

sup ||</j(t)||oo < 2||(/5 o||oo 

0<t<T 

+ C'n( sup ||/+(s)||i+ sup ||<p(s)l|i)l|Vu||So sup ||/+(s 


(4.13) 


(4.14) 




The constant is proportional to ln(T+l) in dimension n = 2 and independent 
of time for u > 3. 

Proof. Fix r > e > 0 and consider t G [£, T]. We shall use the Duhamel formula 
for the heat equation (14.21) with the source term 

/ = /-V-(^Vu). 
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We split the corresponding expression as (^ = concur here that 

Ii is the term corresponding to the initial data, I 2 is the term corresponding to 
the source /, and I 3 corresponds to the convection term —V • {(pVv). For Ii, 
we immediately obtain 

0<li{t,x)=[ K{t,x - y)^o{y)dy <\\(po\\oo ■ (4.15) 


For / 2 , we split the time integral as follows 

rt p (‘t—e 

/O JR'^ 


l 2 {t,x)= / / K{t-s,x-y)f{s,y)dyds= / / ...dyds+ / ...dyds. 

Jo JR" Jo JR" Jt-s Jr" 


■ = Il 


:= 7 ? 


A direct computation shows that 
/ ln(t/£) if n = 2 \ 

\l2\<Cn[ X sup ||/+(s)||l =: C'„$(£,t) sup ||/(s)||l, 

\ £ 2 “ if n > 3 / 0<s<t-e 0<s<t-e 

while 


1 ^ 2 1 < [ if K{t-s,x-y)dy]ds sup ||/+(s)||oo 

Jt—e y JR'^ j t—s<s<t 

< f Ids sup |l/(s)||oo=e sup |l/+(s)||oo. 

J t—€ t—e<s<t t—e<s<t 


As a consequence, we get 


l 2 (t,x)<Cn 4>(£,t) sup ||/+(s)||i+£ SUp ||/+(s)||oo • (4.16) 

\ 0<s<T t-e<s<T j 


Additionally, integration by parts implies that 

rt 


{t,x) = - / K{t-s,x-y)Vy ■ {ip{s,y)Vv{s,y)) dyds 

Jo Jr'^ 

n VyKit- s,x-y) ■Vv{s,y)(p{s,y) dyds 

.n 


/O JR'^ 
pt—e 


/ / ...dyds+ / 

lo JR" Jt-s. 


...dyds =: I 3 + . 
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Similar arguments lead to 


\l3(t^x)\ = 

< 

and 

< 


2i4Tr)^/^ Jo (t-s)(-+i)/2 J 

e~^/^||V^^||oo supo< s<T lly^(5)lli 

y2(47r)"/2(n-1) e(’^-i)/2 


_\x-^ (x — y) 
4 (t-s) ^ 


Vit-s) 


■ Vu (fi dy ds 


(4.17) 


2(47r)"/2 (t _ s)(n+l)/2 

2||e-l"l'x||i||Vu||oo ^ , 

-- y/s sup ||v5(s) 

7r"''2 t-e<s<T 


1 f l^-yd (x — v) 

e — , ■'Vv If dyds 


Vit - s) 


Gathering estimates (14.171) and (I4.18P it follows that 


Thus, combining (|4.15l) . (|4.16|) and (|4.19l) we obtain 
sup ||(^(t)||oo = sup \\h{t) + hit) + hit)\\^ 

e<t<T e<t<T 


(4.18) 


hit,x) < SUp^_^|l(/5(s)||oo^ ■ (4.19) 


^Ilv^olloo+Gn $(e,T) sup ||/+(s)||i+£ sup ||/+(s)||oo 

\ 0<S<T 0<S<T J 

, ^ „ /^supo<,<tI|7>(s)I|i , r- ii ^ mi ^ 

+ G„||Vu||ool --+ Ve^sup^||(/5(s)||oo I- 

Now assume t G [0,e]. A simpler, yet analog, procedure shows that 

sup ||(^(t)||oo = sup Whit) + hit) + l3it)\\oo 
OKtKe 0<t<£ 

< ||</5o||oo +e sup ||/+(s)||oo + ||Vu||oo \/e sup ||(p(s)||oo- 

0<S<T 0<S<T 

Noticing that 

sup ||(/?(t)||oo < sup ||(/?(i)l|oo + sup |l(/3(t)||oo , 

0<t<T 0<t<E e<t<T 


(4.20) 


(4.21) 
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we can add inequalities (|4.20p and (|4.2ip to obtain 


sup ||oo ^ 2|[cpo||oo “t“ Cn 

0<t<T 


^{s,T) sup ||/+(s)||i+e sup ||/+(s)||c 

, 0<s<T 0<s<T 


+ C'nII Vu||oo 


/ SUPo<^<t||‘/?(s)||i 

1 e(»-i)/2 


+ ^/e sup 
o<s<r 



(4.22) 


Let us choose in all dimensions 


e = min {l, 6 ^, T } 


1 + sup 

\ 0<s<T 


\\f+is)\\c 


2 

n + 1 


with 1/5 := 2C'„|| Vu||oo, to discover that 
sup |lv?(t)||oo < 2||(^o||oo 

0<t<T 


+ C'n( sup ||/+(s)||i+ sup ||v?(s)l|l)l|Vu||So +( 1^2 sup ||/+(s)||<;i,+ ' . 
^0<s<T 0<s<T ' 0<s<T 

The constant inherits the dependence of T from $(e,T). Thus, for any 
T > 1 , we get 

Ci(ln(T + l)) ifn = 2, 

( Independent of T if n > 3. 


This proves the lemma. 


□ 


4.2 Analysis of the system fISPDD in dimension n = 2 

4.2.1 Prom to integrability 

Proposition 4.4. Fix a time T > 0 and let (u, w) be a classical nonnegative 
solution of the SPD system (ISPDl) in the interval [0,T]. Then, for any 7 > 1 
there exists an explicitly computable exponent (3 := /3(7) such that 

h(t)ll7 + IkWIU < (1 + , te (o,r], 

with constant C depending additionally o?r ()Hl) -( |HTl ) and || Vpo| 1 l 2 ( 7 +i) but in¬ 
dependent ofT. If additionally (uo, Wo) G x for some 7 > \/2, it follows 
that 

sup (\\u{t)\\-^ + ||u;(t)||^+J < C{mo,'y, ||uo||^, ||wo||.y+) , 

0<t<T '' ' 

where, as before, the constant C depends additionally ort (IH1) -( |HT| ) and || Vpo||l 2 ( 7 +i) 
but it is independent of T. 
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Proof. For the population u, we obtain the following estimate 


— / V? dx + 4, 

dc Jan 


7-1 
7 Jsx 




<( 7 — 1 ) / Vp • V(u"^) dec + 7 / NwvP ^ dec . 

JB" JR" 

With Cauchy-Schwarz and Young inequalities, we get 

(7 — 1) [ Vp • V(u^) dx = 2(7 — 1) f Vp • dec 

J]R'^ 

< 7 ( 7 — 1 ) f |Vp|^u'’'dx + ^ f |V(u'’'/^)|^ dx , 

jR" 7 7 r" 


so that 


4 


vP dx+3 


7-1 


|Vu^/^pdx 


<7 / Nwu^ ^ dx + 7(7 — 1 ) / |Vp|^u'’'dx. 

JR" JK" 


We make use of the Holder inequality with conjugate exponents 7 and 7 ' = 
7/(7 — 1 ) in the first integral of the right hand side and with the pair 7 + 1 , 
(7 + l )/7 in the second. Combined with the convexity inequality ab < dP/p + 
bP /p', we find, for any e > 0 , 


^ / uT'dx+3^^ —- 
at jRn 7 


|Vu'>'/2pdx 


< 7l|lV||oo [ dx + 7 ||iV||oo [ u"^ dx + 

+ [ u'^+^dx. 

7r" Jk" 

In order to control Vp we use estimate (14.31) in Lemma 0 with n = 2, <7 = 
2(7 + 1), 0 = and f = Pw, to conclude that 

[ |Vp(t,x)|^^^’^^^dx 

Jr-^ 

< l|Vpo|| 2 (^+i)+C'(mo)(l + ^) sup ( / |u;(s,x)|^dx) 

|<s<t 

holds for any 7 > v^- The constant C{mo) depends on IHl) and (|H+t . This 
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leads to 


_d 

dt 



u^{t, x) dx + 3-^^- 

7 


a:)p dx 


< 7ll-^lloo / w'^{t,x)dx + 
Jsx- 


7ll^l|oo / V?(t,x) da:+ (l + sup (f w'*{s,x)dxY ^ 

7b" £ ^ ^ ^ |<s<t ^ 

+ 7^e^ x) dx + ^:^|| Vpo|| 2 (;^;J;J|. 

(4.23) 

We also have the following estimate for w 

4- f w“(t,x)dx + 4^^ f \Vw^/^{t,x)\‘^dx (4.24) 

df a Jrti 

< ||c||oo / u“(<,a::)dx + a(||c||oo + ||Aw||oo) / w°‘{t,x)dx. 

In other to control the right side with the left side we are forced to select 7 < a < 
7 + 1. We shall also use Gagliardo-Nirenberg-Sobolev’s inequality: in the whole 
space inequality (13.41) reduces to J dx < C J ^ dx / V(^“/^) dx. Thus, 
adding inequalities (14.231) and (|4.24l) we obtain, after similar computations to 
those of the FPD system, the following estimate which holds for any t G (0, T] 

-^Z(t)+ <7(7710)^“'(t) < C(l + ||Vpo|l 2 |;^+i])+C(mo)(l + ^) ^sup II^II^T'' 

2 <s<t 

< C(1 + + C{mo){l + ^sup ||w||;(“' 

2 —^ 

< (7(1 + \\ypo\\l4-i]) + C'(’™o)(l + ^sup Z-“'(t). 

2 —^ 

where 

Z{t) := f M^(<)dx+ f w°'{t)dx. 

7r"- JRn 

Now apply the comparison Lemma lA.21 to obtain for any 7 > y/2 

Z{t) <C{mo,-f)(l +, 0 <t<T. (4.25) 

V ^ 7 ^ 7 / 

The constant C depends additionally on(IHl)-( |H+D and ||Vpo||L 2 (-y+i) but it is 
independent of T > 0. The case 7 G (1, \/2] follows by Lebesgue’s interpolation 
between estimate (j4.25|) and the mass conservation. 
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Finally, uniform propagation of the L'* and norms of u and w follows 
using previous computations with the estimate (14.41) instead of estimate (I4.3|) 
which give us the bound 

Az{t) + Cimo)Z‘^\t)<C{l + \\Vpo\\l^r^ll])+C{mo) sup 

w-r ; 0<s<T 

valid for any t € [0, T] with T > 1 and \/2 < 7 <a <7 + l. The comparison 
Lemma lA.ll gives 

sup Z(s) < 

0<s<T 

1 

max|z( 0 ), (^(7(1 + || + ^(mo) ^sup^Z°“'(s)) “ |, 

which implies, since ^ < 1, that supQ<g<y Z(s) is finite and uniform in T. □ 

4.2.2 Prom L'* to L°° integrability 

Proposition 4.5. Let T > 0. Consider initial data (^Uo, Wo, ypo) G x x 
L® and, let {u, w) be a classical nonnegative solution of the system (ISPDl) in 
the interval [0,T]. Then, the following L^-estimate holds 

||u(t)||oo + Ik(t)||oo < , fe(0,T], (4.26) 

where the constant C depends on (0, dig and the initial data, but, it is 
independent of T > 0. In particular, for any 7 G [1, 00 ] it follows that 

||u(t)||^ + ||w(t)||^ < , te(0,T]. (4.27) 

Furthermore, estimate (|4.26|) can be upgraded by adding the dependence on the 
L°° -norms of the initial data in the constant, 

sup ||t/;(5)||oo + sup ||w(0lloo < (^(mo, ll'u^olloo, ||uo||oo) , (4.28) 

0<s<T 0<s<T 

The constant depend on (IHl, dig and ||Vpo||oo; but, it is independent of the 
time. 

Proof. The proof of estimate (14.261) is a direct consequence of De Giorgi level set 
technique presented in the proof of Proposition 13.31 and Proposition 14.41 The 
integrability of the initial data assures that the delicate nonlinear term u Vp 
will remain uniformly bounded as required in these arguments. 

Interestingly, the uniform bound (14.281) is not straightforward to prove. The 
main reason is the lack of knowledge about the regularity of Ap, which was a 
consequence of elliptic regularity for the FPD system. This is where we make 
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use of Lemma 14.31 Indeed, let us apply Lemma 14.31 to the population u, that is, 
if = u, f+ = w N, and recalling that conservation of mass implies 

sup ||w(s)||i+ sup 11 ^( 5)111 < 2mo . 

0<S<T 0<S<T 

Thus, we get 

sup ||'u(s)||oo ^ ^IIUolloo 
0<s<T 

^ (4.29) 

+ Cl sup ||Vp(s)||((^+ sup |lw(s)||<i+'. 

0<s<T 0<s<T 

Now, apply Lemma |4.11 with q = 00 and 0 = 2, to the pheromone equation to 
obtain 

sup ||Vp(s)||oo < 2||Vpo||oo + sup ||■lc(s)||<)(^ . (4.30) 

0<s<T 0<s<T 

Noticing that c < Co, we may apply Lemma 14.31 to the population w as well to 
conclude that 

sup ||u;(s)||oo < 2||ri;o||oo+ C'^||Vu||((„ sup ||u(s)||^ . (4.31) 

0<s<T 0<s<T 

In addition to the mass, the constants depend on the L°°-norms of the param¬ 
eters and grow logarithmically in time. Plugging successively Ij4.31|l in Ij4.30p 
and the result in (j4.29ll the estimate for the population u reduces to 

sup \\u(t)\\oo <Co + Cl sup ||u(s)||oo"^^ + C 2 sup ||u(s)||(4.32) 

0<t<T 0<S<T 0<S<T 

where Co depends on the initial data and ||Vu||oo- Since | < 1 the 

result follows for short T > 0 with constants growing at most ln(T). For large 
T > 0 the result follows after simple interpolation with estimate (I4.26P . □ 

Remark 4.6. The maximal exponent in the right side of (14.321) > 1 for 

n > 3. This suggest, for these dimensions, a finite in time Dirac concentration 
of mass similar to that of the Keller-Segel model in two or more dimensions. 
Similarly, such concentration could be avoided by smallness conditions on the 
initial data and model parameters. 


A A useful comparison lemma 


Lemma A.l (ODE comparison). Assume Y and X are absolutely continuous 
functions in [0, T] and such that 


Y'{t)+aY°‘{t)>h + 5 + c(l + sup r“°(s) 

C ^ T<S<t 

X'{t) + a X°‘{t) < b + c(l +—\ sup X““(s), 

V t')'/ T-<s<t 


(A.l) 
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with 6 > 0, c > 0, a > 0, (5 > 0, a > ao > 0, 7 > 0 and t > t > 0. If 
y(0) > X(0) then Y > X in [0,T]. In particular, 1/7 = 0 

sup X(i) < max{X(0), C} , (A.2) 

tG[0,T] 

where the constant C > 0 depends on all parameters but t, S and T. 

Proof. Define 

to := sup : D(s) > X{s), s S [0,t]} . 

Note that to > 0 since 1^(0) > ^(0). Let us argue by contradiction assuming 
that there exists ti G (0,T] such that X{ti) > Y{ti). Since X and Y are 
continuous to G (0,T) and X(to) = Y(to). Also, there exists interval I = 
{to, to + s) such that X{t) >Y(t) for any t G I. Thus, the fundamental theorem 
of calculus implies 

^ X'{s) ds = X{t) - X{to) > Y{t) - Y{to) = Yfs) ds , tGl. 

Therefore, there exists t.^ G I sufficiently close to to such that: (1) A" and Y are 
differentiable at with > Y'{p), and (2) c(l + ^) (supT-<g<(^ Y°'° (s) — 

supT<j,<t^ Ar“‘’(s)) > —S. Thus, using (lA.ll) 

0 > Y'(U) - X'(U) > a(X^(U) - F“(t*)) 

+ (5 + c(l + i')( sup y“°(s)— sup X“°(s)')> 0 . 

This contradicts the existence of ti. Finally, the estimate (IA.2I) follows by taking 
Y := Ys as the constant function in [0, T] given by max{Ar(0) + 5, Cs} where 


Cs = -{h + 5+ £“/“« = 


a 

4c 


Then, X{t) < Ys for any 5 > 0. The result follows by sending 5 —0. 


□ 


Corollary A.2. Assume X be an absolutely continuous function in [0,T] such 
that 

X' + aX°‘ <b + c(l + ^') sup X“°(s), 

with b > 0, c > 0, a > 0, a > ao > 0 and 7 > 0. Then, 

X{t)<c[l + ^'^, /3 = max{^,^}. 

The constant C > 0 depends on all parameters hut it is independent of T. 
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Proof. First note that the function 


y{t) = c(i +1) 

satisfies Y'{t) + aV^it) > & + <5 + (l + ^) sup y“°(s) for C sufficiently 
large depending on all parameters. Second, since X is bounded on [0,T], there 
exists sufficiently small to > 0 such that y(s) > supjgfg t] * G (0,to). 

Applying Lemma [A. II in the interval [s,T] it follows that Y > X in [s,T]. The 
result follows since s can be taken arbitrarily small. □ 
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